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Until  recently,  the  capability  to  accomplish  a two-dimension»l  axisymmetric  nozzle  thermal 
analysis  wnieh  considers  in-depth '< charring  has  not  been  available/€o  rfchty  rocket  nozzle  analysis 
^ community- except  by  using  "pseudo1’-' material  properties  to  account  for  in-depth  decomposition  and 
using  a standard  2-D  axisymmetric  conduction  code  with  surface  ablation. 

To  alleviate  this  problem  a two-dimensional  axisymmetric  heat  transfer  and  ablation  computer 
code  has  been  developed  which  accounts  for  in-depth  decomposition.  The  code  was  developed  by 
extensively  modifying  the  Acurex  ASTHMA  II  code  to,  account  for  in-depth  charring,  grid  generation 
and  material  element  conductivity  based  upon  the  local  element  coordinate  system.  The  code  was 
also  changed  to  handle  as  many  as  2,000  elements.  C'cclf®  <’(<>  ,•  V s n . ,'t"  ► <■' ••  ■ /■.■*■  t -}  o-,‘ 

The  development  of  the  code  was  based  upon  modeling  the  material  response  to  a thermal 
environment  in  the  same  way  as  the  one-dimensional  Acurex  CMA  program.  Example  problems  are  shown 
to  compare  the  charring  version  of  ASTHMA  (ASCHAR)  to  typical  CMA  models  showing  the  resulting 
ablation  and  thermal  profiles.  3 , 

Example' problems  showing  the  results  of  the  ASCHAR  analysis  of  a complete  nozzle  will  be*  A 1 * 
given  showing  two-dimensional  thermal  effects.  _ 

K 

INTRODUCTION 


Heretofore,  the  rocket  motor  design  community  could  not  adequately  analyze  rocket  nozzles 
because  of  the  lack  of  a 2-D  axisymmetric  heat  transfer  code  which  models  in-depth  material 
decomposition  and  the  subsequent  energy  transport  due  to  the  pyrolysis  gases  percolating  through 
the  char  layer  to  the  surface.  In  general,  most  of  the  rocket  motor  community  used  Aerotherm's 
(Acurex)  1-D  "Charring  and  Ablation  Code"  (CMA)  to  model  charring.  2-D  codes  were  available, 
but  most  of  the  codes  which  were  being  used  did  not  have  charring  or  material  decomposition  capab- 
ilities and  as  a result  most  of  the  solid  rocket  motor  manufacturers  used  these  noncharring  codes 
with  pseudo  material  properties  to  account  for  charring  when  2-D  analyses  were  required. 

The  disadvantages  of  using  a noncharring  2-D  code  with  pseudo  properties  are: 

1.  It  is  difficult  to  adequately  model  the  properties  in  all  temperature  ranges. 

2.  Different  heating  rates  cannot  be  accounted  for,  especially  if  at  some  time  the  cooling 
effects  are  modeled  such  as  in  post  test  cooling  and  in  pulse  motors  when  there  is  a 
nonheated  coast  period. 

Because  of  the  lack  of  an  adequate  tool  to  determine  the  2-D  charring  effects  in  rocket  motor 
nozzles,  Morton  Thiokol  began  to  develop  such  a code  in  the  mid-1970s.  It  was  determined  at  that 
time  that  the  most  efficient  means  of  developing  such  a code  would  be  to  modify  the  original 
Aerotherm  "Axisymmetric  Transient  Heating  and  Ablation"  (ASTHMA)  code.  The  main  reasons  for  using 
this  code  as  a basis  were: 

1.  The  code  was  already  recognized  in  the  industry  for  nozzle  thermal  analysis  expecially 
for  nozzles  with  graphite  and  carbon/carbon  throats  and  ITEs. 

2.  Heterogenous  surface  chemical  reactions  were  modeled  for  surface  ablation  calculations. 

Even  though  the  ASTHMA  code  was  a good  basis  to  start  the  development  of  a 2-D  charring  code, 
extensive  modification  was  required  because  of  the  following  limitations. 

1.  Math  modeling  and  numerical  development  for  charring  was  necessary. 

2.  ASTHMA  required  that  all  elements  (nodes)  in  a material  had  to  be  oriented  in  the  same 
direction.  (No  possibility  of  contouring  the  elements  with  material  boundaries.) 
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3.  Original  version  of  ASTHMA  was  for  200  elements  maximum  and  as  a result,  a complete  nozzle 
could  not  be  analyzed. 

4.  There  were  limitations  on  applying  different  heating  rates  to  many  surface  elements  of  a 
complete  nozzle  because  of  limitations  on  the  number  of  boundary  condition  tables. 

BACKGROUND  AND  THEORY 

T«o  types  of  thermal  analyses  are  generally  done  for  estimates  of  temperatures  and  material 
ablation  in  rocket  nozzles.  The  first  is  a one-dimensional  analysis  which  is  generally  mode  for 
preliminary  estimates  of  nozzle  temperatures  and  also  the  prediction  of  the  surface  and  indepth 
ablation.  The  analysis  is  completed  by  analyzing  the  nozzle  at  several  critical  locations  such  as 
the  throat,  nose,  and  exit  cone. 

A two-dimensional  (2-D)  analysis  is  generally  made  to  determine  the  effects  of  irregular 
geometry  and  to  obtain  a detailed  map  of  the  temperatures  in  the  nozzle  so  that  a detailed  struc- 
tural analysis  can  be  made.  The  scope  of  this  paper  only  addresses  the  2-D  analysis  and  as  a 
result  the  background  and  theory  necessary  for  this  type  of  analysis  will  be  discussed. 

The  basis  for  the  2-D  charring  ASCHAR  code  is  essentially  the  same  as  Aerotherm's  (Acurex) 
"Charring  Material  Thermal  Response  and  Ablation  Program"  Version  3 (CMA)  except  the  ASCHAR  is  a 
2-D  axisymmetric  code.  Figure  1 illustrates  the  general  physical  problem  treated  by  the  ASCHAR 
code. 
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Figure  1.  General  Physical  Froblem 


The  surface  energy  and  mass  balance  considered  by  the  ASCHAR  code  is  depicted  in  Fig.  2. 


Figure  2.  Surface  Energy  and  Mass  Balance 
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This  energy  balance  can  by  described  mathematically  as: 


PeVi^r  - hew>  + PeueCM  (Z!w)hiw)  ' B'hw  + mchc  + %h 


8 8, 


’sen 


’chem 


w’rad 


- Fae  T 


w w 


’cond 


= 0. 


(1) 


-J  ^ 


’rad 

in 


qrad 


out 


where: 


Pe  = local  gas  density  at  the  edge  conditions 


u = local  gas  velocity  at  the  edge  conditions 


Cjj  = Stanton  number  for  heat  transfer 
= recovery  enthalpy 

hgw  = enthalpy  of  the  edge  gas  at  the  wall  temperature 
Cjj  = Stanton  number  for  mass  transfer 

Z*e  = fraction  of  specie  i at  the  edge  due  to  chemical  equilibrium  and  diffusion 

Z*w  = fraction  of  specie  i at  the  wall  due  to  chemical  equilibrium  and  diffusion 
T 

h^  w = enthalpy  of  specie  i at  the  wall  temperature 
B'  = dimensionless  blowing  rate 
= (*,  + Ac)/PeueCM 

8 


m_  = mass  flow  rate  of  pyrolysis  gas  out  from  surface 


mc  = mass  flow  rate  of  char  from  surface 


hw  = enthalpy  of  the  wall  gases 


CRA&I 
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h - enthalpy  of  char  at  wall  temperature 


hg  = enthalpy  of  pyrolysis  at  wall  temperature 


■oution  I 


<*w  = surface  absorptivity 


qra(j  = incident  radiation  flux 


F = view  factor 


emissivity  of  wall 


Tw  = wall  temperature 


qcond  = conducted  into  wall 


The  first  term  of  equation  (1)  represents  the  sensible  energy  convected  to  the  surface,  the 
qchemterni  represents  a number  of  chemical  energy  fluxes  at  the  surface  and  the  remainder 
represents  the  radiation  and  conduction  of  energy  to  and  from  the  wall. 


The  ASCHAR  computer  program  solves  the  surface  energy  balance  equation  in  much  the  same  way 
as  the  original  (ASTHMA  II)  code.  However,  in  the  indepth  en  -gy  balance,  the  original  version  of 
the  program  does  not  have  the  capability  of  internal  decompose .ion  and  the  reference  coordinate 
system  is  not  at  the  wall  surface.  As  a result  no  terms  are  associated  with  the  pyrolysis  gas 
leaving  the  surface  (a  h ) or  energy  transfer  due  to  coordinate  system  change. 


253 


85363-1.4 


ASCHAR  has  this  capability  of  indepth  charting,  and  as  a result  the  energy  associated  with  gas 
injected  into  the  surface  control  volume  is  necessary  and  a term  was  added  to  account  for  this 
energy  Likewise,  additional  terms  were  added  to  the  indepth  conduction  solution  to  account  for 
the  energy  associated  with  decomposition  and  convective  energy  associated  with  pyrolysis  gas 
moving  through  the  material  to  the  surface. 

The  equation  for  the  indepth  conduction  solution  is  an  explicit  finite  difference  type  as 
opposed  to  the  implicit  solution  employed  by  CMA.  The  equation  is: 

T ^m+1  ,n,6  ^tn.n+l  ,6  ^m-1  ,n,8  . ^m,n-l  ,8  _ 

“ R n Q R .C.  R,  „ a R , A Q ' m,n,6 
u,n,B,8  m,n,A,0  m-l,n,B,0  m,n-l,A,8  y 

i + i — + i — + i 
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pyr  m.n,0  g m-l,n,0 

dH 


- M 
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l dT 

m-l,n,0 

A0 
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where:  T = temperature  of  element  m,  n at  the  new  time  0’ 

m,n,o 

R = thermal  resistance 
A = side  A of  element 
B side  B of  element 
0 = old  time 

M a = generated  gas  rate  of  element  m,n 

pyr  m,n,o 

Mm  n 0 = tota^  8as  rate  in  element  m,n  (generated  plus  gas  from  element  m-l,n) 
Hg  = enthalpy  of  the  pyrolysis  gas 

H = the  equivalent  enthalpy  of  the  material  that  changes  from  solid  to  gas 
defined  as  (pp  hp  - pc  hc)/(pp  - Pc) 

pp  = density  of  virgin  plastic 

Pc  = density  of  char 

hp  = enthalpy  of  virgin  plastic 

hc  = enthalpy  of  char 


dH 
dT 

A0  = time  increment 


— =£  = slope  of  enthalpy  vs  temperature  curve  for  pyrolysis  gas 


C = thermal  capacitance  pVCp 

Note  that  the  energy  associated  with  the  pyrolysis  gas  pickup  is: 

m dH  T . M dH 

m-l,h,0  dT  m-l,n,0  m-l,n,0  m,n,0  dT  m,n,0  m,n,8 

and  the  term  for  the  material  decomposition  is: 


Mpyr  m,n,6  ^Hg  " 
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Ap  = -B. 
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These  terns  are  not  included  m the  original  ASTHHA  program. 


Another  difference  between  ASCHAR  and  ASTHMA  is  in  the  way  internal  thermal  conductivity  is 
treated.  The  basic  heat  conduction  model  employed  in  both  programs  is  based  upon  one-dimensional 
flow  between  adjacent  elements.  Consequently,  anisotropic  heat  transfer  cannot  be  modeled  exactly 
except  in  the  limited  case  in  which  the  axes  of  every  element  coincide  with  the  material  properties 
axes.  However,  ASCHAR  provides  an  approximate  solution  when  these  axes  do  not  coincide.  Figure  3 
shows  how  this  is  accomplished.  As  shown,  k^  and  k„  are  the  conductivities  along  the  material 
major  and  minor  axes,  and  0-c  and  0^^  are  the  orientations  of  the  element  axes  relative  to  the 
material.  ASCHAR  approximates  the  conductivity  in  the  AC  and  BD  directions  from  the  following: 


kAC  = kj  cos2  eAC  + k2  sin2  eAC 


(A) 


kBD  = kl  COs2  0BD  + k2  Si“2  0BD 


Figure  3.  Material  Conductivities  Solution 


Boundary  conditions  and  surface  thermochemistry  used  in  the  ASCHAR  program  are  calculated  by 
such  programs  as  ACE  or  GASKET  and  MEIT.  See  References  4,  7 and  8. 


ASCHAR  INPUT 


Basically  there  are  four  major  input  efforts  required  to  run  the  ASCHAR  program.  They  are 
the  noz2le  geometry,  the  element  (node)  data,  the  thermochemistry  data  and  the  boundary 
conditions.  Most  of  the  input  is  straight  forward  and  the  input  document  is  all  that  is  required. 


2SS 


Edge  and  surface  thermochemistry  tables  are  required  by  ASCHAP  Thp  ACE  computer  code  is 
generally  used  to  generate  this  input,  however  other  programs  such  as  GASKET,  or  hand  generated 
tables  are  also  acceptable.  Several  thermochemistry  tables  at  different  pressures  are  selected 
for  different  locations  around  the  nozzle  as  shown  in  Fig  k which  describes  one  of  the  example 
problems  that  was  analyzed. 


The  input  for  the  element  data  or  "node"  data  requires  some  discussion  about  the  method  of 
computing  the  multiplying  factors  for  the  heat  transfer  coefficient. 

Because  only  a limited  number  of  boundary  condition  ("function  of  time")  tables  are  allowed, 
the  program  was  modified  to  accept  multiplying  factors  on  the  heat  transfer  coefficient  so  that 
the  heat  transfer  coefficient  could  be  changed  from  node  to  node  in  the  region  that  applies  to  a 
particular  "function  of  time"  table.  As  shown  in  Fig.  5,  there  were  eleven  "functions  of  time" 
tables  assigned  to  different  regions  of  the  nozzle  for  the  carbon  phenolic  exit  cone  nozzle. 


The  "function  of  time"  tables  used  in  ASCHAR  are  required  for  the  recovery  enthalpy, 
radiation  flux  and  heat  transfer  coefficients.  A discussion  of  how  these  variables  are  obtained 
is  presented  in  the  following  paragraphs. 

The  recovery  enthalpy  is  taken  directly  from  the  plot  for  the  recovery  enthalpy  of  MEIT  output 
(Fig.  6). 


Figure  6.  Carbon  Phenolic  Nozzle  Recovery  Enthalpy  Distribution 
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The  radiation  input  for  the  "function  of  time"  tables  is  computed  based  upon  the  equation: 


arr ' 


where: 


Qt  = incident  radiation  heat  flux  (Btu/ft  -sec) 


o = Stefan-Boltzmaan  constant  (4.7611  x lO-13  Btu/ft2-sec-°R^) 


Tt  = local  average  static  temperature  of  combustion  products  (°R) 
e = emissivity  of  local  combustion  products,  1-exp  [-0.808  (%A1 / 16)  D] 
%A1  = percent  aluminum  in  propellant  composition 


p = local  density  of  combustion  products  including  condensed  species,  (Ibm/ft  ) 
D = local  nozzle  diameter  (in.) 


The  radiation  flux  is  not  generally  made  a function  of  time  because  the  flame  temperature  and 
resulting  edge  temperatures  do  not  change  appreciably  with  pressure;  therefore,  the  value  is  used 
as  a time  average  value.  This  value  is  used  for  the  entire  pressure  region;  and,  therefore,  a 
surface  node  should  be  selected  near  the  center  of  the  region  so  that  the  value  will  represent  an 
average  value  for  the  entire  region  (see  Fig.  5). 


Figure  7.  Comparison  o:  ASCHAR  with  Cf*A  in  a Carbon  Phenolic  Nozzle  Inlet 


The  heat  transfer  coefficient  is  made  a function  of  time  for  a node  in  a particular  region  of 
the  following  relationship: 


p U C.  , H (6)  = p U C,  , N 
Ke  e h’  'e  e h’ 


P(6) 


0.8 


where:  PeUeCh’  ^ ^ = heat  tran;:*e,:  coefficient  of  a node  in  a region  at  time  0 


Pe'eSi’  ^ = avera8e  heat  transfer  coefficient  cf  a node  in  a region 
P(8)  = pressure  at  time  6 

p = average  chamber  pressure 


The  ratio  of  mass  transfer  to  heat  transfer  (C  / C is  supposed  to  represent  the  ratio  of 
Stanton  number  for  mass  and  heat  transfer.  This  vaTue'is  approximated  by: 


> =Le2'3=0.7 
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However,  for  most  materials  except  carbon  phenolic,  this  factor  can  be  used  as  a correlating  factor 
to  correct  the  "theoretical"  ablation  to  match  observed  rocket  motor  firings.  A value  of  C^/C^  = 
0.3  has  been  reasonably  successful  when  using  ACE  thermochemistry  for  carbon/carbons  similar  to 
that  used  in  the  example  nozzles  (figs. 8 and  9).  'fhis  value  is  shown  in  the  lead  card  of  the  sur- 
face thermochemistry  data  which  is  generated  by  ACE. 

ASCHAR  VERIFICATION  AND  NOZZLE  THERMAL  ANALYSES  RESULTS 


To  verify  that  the  ASCHAR  code  is  a reliable  tool  for  nozzle  thermal  analysis,  the  code  was 
checked  by  doing  simple  1-D  analyses  using  charring  ablators  and  comparing  the  resulting  thermal 
profiles  .'ith  the  1-D  CMA  results.  Figure  7 shows  the  results  of  a typical  analysis  and  compares 
the  2-D  ASCHAR  with  CMA  using  the  same  material  properties  and  boundary  conditions. 

To  verify  that  the  code  can  actually  perform  complex  2-D  axisynxnetric  nozzle  thermal 
analysis,  the  results  of  the  thermal  profiles  and  surface  ablation  of  three  different  types  of 
nozzles  as  computed  by  the  ASCHAR  code  are  shown  in  Figs.  8 through  10. 


Figure  8.  Analysis  Results  of  Strategic  Nozzle  with  C carbon  Phenolic  Exit  Cone 
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Figure  9.  Analysis  Results  of  SRM  Nozzle 
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